A note on limit shapes of minimal difference partitions 
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Abstract 

We provide a variational derivation of the limit shape of minimal differ- 
ence partitions and discuss the link with exclusion statistics. 

This paper is dedicated to Professor Leonid Pastur for his 70th anniversary. 

A partition of a natural integer [1] is a decomposition of as a sum of 
a nonincreasing sequence of positive integers {hj}, i.e., E = Ylj^j such that 
hj > /ij+i, for j = 1, 2 . . .. For example, 4 can be partitioned in 5 ways: 4, 3 + 1, 
2 + 2, 2 + 1 + 1, and 1 + 1 + 1 + 1. Partitions can be graphically represented 
by Young diagrams (also called Ferrers diagrams) where hj corresponds to the 
height of the j-th column. The {/ijj's are called the parts or the summands of the 
partition. One can put several constraints on such partitions. For example, one 
can take the number of columns to be fixed or put restrictions on the heights. 
In this paper we focus on a particular constrained partition problem called the 
minimal difference p partitions (MDP-p). The MDP-p problem is defined by 
restricting the height difference between two neighboring columns, hj — /ij+i > p- 
For instance the only allowed partitions of 4 with p = \ are 4 and 3 + 1. A typical 
Young diagram for MDP-p problem is shown in figure [H Consider now the set of 
all possible partitions of E satisfying E = hj and hj — hj+i > p. Since this is 
a finite set, one can put a uniform probability measure on it, which means that 
all partitions are equiprobable. Then, a natural question is: what is the typical 
shape of a Young diagram when E — > cxd? 

In the physics literature this problem was first raised by Temperley, who was 
interested in determining the equilibrium profile of a simple cubic crystal grown 
from the corner of three walls at right angles. The two dimensional version of 
the problem — where walls (two) are along the horizontal and the vertical axes 
and E "bricks" (molecules) are packed into the first quadrant one by one such 
that each brick, when it is added, makes two contact along faces — corresponds 
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Figure 1: A typical Young diagram for MDP-p problem. The thick solid border 
shows the height profile or the outer perimeter. Wh is the width of the Young 
diagram at a height h, i.e., Wh is the number of columns whose heights > h. 



to the p = partition problem. Temperley [2] computed the equilibrium profile 
of this two dimensional crystal. In the mathematics literature the investigation 
of the shape of random Young tableaux was started by Vershik and Kerov [3] and 
independently by Logan and Shepp The case of uniform random partitions 
was treated by Vershik and collaborators [5l[6l[7j who obtained the limit shapes for 
the p = and p = 1 cases and also the average deviations from the limit shapes [8] . 
Some of these results were extended by Romik [9] to the MDP-p for p = 2. These 
problems belong to the general framework of asymptotic combinatorics, a subject 
which displays unexpected links with random matrix theory [lOj . In this note we 
compute the limit shapes of MDP-p for all p > by a variational approach and 
mention an interesting link with exclusion statistics. 

We first recapitulate the arguments used in [111 Il2j (see also [13] for a similar 
approach) to compute the limit shapes of the Young diagrams of unrestricted 
partitions (p = 0). Let P = {i,hi) and Q = (j, hj) be two points belonging to the 
outer perimeter of the Young diagram of a given partition. We evaluate the total 
number of subdiagrams which connects these two points. These subdiagrams 
are lattice staircases with the only restriction that each step either goes right or 
downward. The total number of horizontal steps is j — i, the total number of 
vertical steps is hi — hj , and the total number of steps \s j — i + hi — hj . Therefore, 
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the total number of configurations is 

no{P,Q) = no{i,hi;j,hj) 



j-i + hi - h. 



If P and Q are far apart (i.e., a = j — i » 1, 6 = /ij — /ij ^> 1) we may use the 
Stirhng formula which gives 



\n9.Q{P,Q) = -aln- 
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a + b a + b 

where Ti = (ni, n2) = {b, a)/ Va?~+1^ is the unit vector orthogonal to PQ and 
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Heuristically one expects that in the limit £" — > oo,/i — > oo, Wh 
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profile of the Young diagram will be described by a smooth curve y = y{x) where 
y = hj^/E and x = Wh/^/E are the scaling variables. The normal vector can be 
parameterized as 
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In the lattice model, the points P and Q were taken to be far apart. However 
in the new scale (x, y) one now assumes that they are close enough in order to 
ensure that the interface is locally flat. The total number of Young diagrams O 
with a given area E is obtained by adding all such local configuration, i.e. 



n = exp(\/^/ dxA/1 + y'2(2;) 0( 
with the area constraint 
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For large the most dominant contribution to arises from the optimal 
curve y = y{x) which maximizes the integral in ^ with the constraint ([6]). This 
optimal curve describe the limit shape of the Young diagrams. Thus we are led 
to the variational problem of extremizing 
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A / y(x)dx, (7) 



3 



where A is a Lagrange multiplier. This leads to the Euler-Lagrange equation, 
which simplifies to 

We solve this equation with the boundary conditions y(oo) = and y{x ^ 0) — > 
oo. The later condition follows from the fact that y = h/\fE ~ Ini? when 
X = Wh/^fE for large E [_14j. Therefore y(0) diverges in the limit E — > oo. 
The solution gives the equation of the limiting shape as 

2/(x) = -^ln(l-e-^-) with A = ^, (9) 

where A is obtained by using the constraint ([6]). 

The goal of this paper is to extend this derivation to the MDP-p with p > 0. 
This is a priori non-trivial since now one has to take into account the restriction 
on the steps. In the following we will use an exact correspondence between 
MDP-p with p > and a unrestricted partition {p = 0). 

Let {/ij} denote the set of non-zero heights in a given unrestricted partition 
{p = Q) E = Yl!j=i^ji where hj > /ij+i for all j = 1, 2, . . . , — 1. Let us 
now define a new set of heights h'- = hj + p{N — j) for j = 1, 2, . . . , A^. Thus 
h'j — h'-j^^ = hj — hj+i + p for all J = 1, 2, . . . , — 1 and h'j^ = h^ > 0. Since 
hj — /ij+i > 0, the new heights thus satisfy the constraint h'j — h'-_^^ > p for all 
J = 1, 2, . . . , — 1. Since the mapping is one to one, the total number of local 
MDP-p configuration satisfies 

9.p{i,h[;j,h'j) = Qo{i,hi;j,hj). 

Moreover, hi — hj = h'^ — h'- — p{j — i). Therefore using ([1]), 

^ ^ ^'[j-i){l-p) + h',-h'; 
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The fact that the mapping does not preserve the total area does not spoil the 
argument since here we only deal with local MDF-p configurations. The area 
constraint is a global one which is implemented at the end of the calculation 
via a Lagrange multiplier. Following the same steps as before we arrive at the 
variational problem of extremizing 



y{x) dx. 



(10) 

Using the same Euler-Lagrange formalism, finally leads us to the equation of the 
limit shape for p > 0, 

y = -iln(l -e"-^'') -px. (11) 
A 
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Figure 2: Limit shapes for the minimal difference p partitions with p = 0, 1, 2, and 
3, where A(0) = tt/VG, A(1) = tt/V12, A(2) = vr/vTS, and A(3) = 0.752617.... 



The Lagrange multipher A in (jlip can be determined by using condition y > 
and the normahzation J^^™ dx = 1, where Xm is the solution of the equation 
y{xm) = 0. Writing exp(xm,) = y*, it satisfies y* — y*^~^ = 1, and in terms of y* 
one finds 

A^ ^ \\p) = ^-U,{l/y*) - |(lny*)2, (12) 

where Li2(2;) = "^2^=1 is the dilogarithm function. X{p) is a constant which 

depends on the parameter p. For instance for p = 0, 1 and 2, one finds A(0) = 
7r/\/6, A(l) = 7r/-v/T2 and A(2) = vr/VTS in agreement with the earlier known 
results [5l [9]. Figure [2] shows the limit shapes for the MDP-p with p = 0,1, 2, 
and 3. 

Equation (fTTj) implies that the inverse function x{y) = ln(j){y) satisfies 

4>{y)-e-^yct>{yy-P = l (13) 

Amazingly this equation appears in several apparently unrelated contexts. 

1. The generating function S(t) = 1 + X^^Li Sk{q)t^ for the number of con- 
nected clusters Sk{q) of size /c in a q-ary tree satisfies [15] 

S{t) - tS^t) = 1. (14) 

This establishes a formal link between two different combinatorial objects, on one 
hand the q-ary trees and on the other hand the MDP-p problem with p = 1 — q. 
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In graph theory |16j . Sk{q) is known as the generahzed Catalan number, which is 
given by 

If qk 



2. Consider the generating function of MDP-p problem 



Z(x,z) = ^j;Pp(i?,iV) 



E N 



where pp{E,N) is the total number of MDP-p of S in parts. In the limit 
S — > oo the number of such partitions will be controlled by the singularities of 
Z{x, z) near x = 1. By setting x = , one gets for /3 — > [T7j 

\nZ{x,z)^J^ Inyp (^ze~^') de, (15) 

where the function yp{t) is given by the solution of the equation 

yp{t)-tyl-P{t) = l. (16) 

3. In the physics literature (I13p also arises in the context of exclusion statis- 
tics. Exclusion statistics \18\ [19\ [201 Ell US] — a generalization of Bose and Fermi 
statistics — can be defined in the following thermodynamical sense. Let Z(P, z) 
denote the grand partition function of a quantum gas of particles at inverse tem- 
perature P and fugacity z. Such a gas is said to obey exclusion statistics with 
parameter < p < 1, if Z{P, z) can be expressed as an integral representation 

In Z(/3, z) = /5(e) In (z^-^^) de, (17) 

where p{e) denotes a single particle density of states and the function yp(t), 
which encodes fractional statistics is given by the solution of (jlSp . Well known 
microscopic quantum mechanical realizations of exclusion statistics are the Low- 
est Landau Level (LLL) anyon model pl9] and the Calogero model [20], with p{e) 
being, respectively, the LLL density of states and the free one dimensional density 
of states. 

The fact that the same equation appears in all three cases is obviously not 
fortuitous. The link between 2 and 3 follows from the fact that exclusion statis- 
tics has a combinatorial interpretation in terms of minimal difference partitions 
which generalizes the usual combinatorial interpretation of Bose statistics (resp 
Fermi) in terms of partitions without (with) restrictions. Let us briefly recall this 
correspondence. Let be the number of columns of height h = i m. a Young di- 
agram of a given partition of E, then E = rijej can be interpreted as the total 
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energy of a non- interacting quantum gas of bosons where Ci = i for i = 1, 2, . . . , oo 
represent equidistant single particle energy levels and = 0,1, 2,... ,00 repre- 
sents the occupation number of the i-th. level. If one now puts the restriction 
that hj > /ij+i (e.g. allowed partitions of 4 are: 4 and 3 + 1), then the restricted 
partition problem corresponds to a non- interacting quantum gas of fermions, for 
which rii = 0, 1. If in addition, one restricts the number of summands to be A^, 
then clearly = ^ - rii represents the total number of particles. For example, if 
E = 4 and N = 2, the allowed partitions are 3 + 1 and 2 + 2 in the unrestricted 
problem, whereas the only allowed restricted partition is 3 + 1. The number 
p{E, N) of ways of partitioning E into parts is simply the micro-canonical 
partition function of a gas of quantum particles with total energy E and total 
number of particles A^: 



For both unrestricted and restricted partitions, one can readily check that the 
grand partition function Z{e~^,z) = ^^^"^^/^(^i limit P — > 

0, is nothing but the one in ()15p . with p = and p = 1 respectively. 

For a quantum gas obeying exclusion statistics with parameter p it is a priori 
not obvious how to provide a combinatorial interpretation since the underlying 
physical models with exclusion statistics describe interacting models. However in 
some specific cases , such as the Calogero model one can show that the spectrum 
can be parameterized as a free spectrum with some restrictions on the quantum 
numbers which reflect the fact that the Pauli principle is replaced by a stronger 
exclusion principle [23\ [23] . This exclusion is enforced at the level of the Young 
diagrams by the constraint hj — /ij+i > p- The link between 1 and 3 expresses 
this correspondence in terms of counting of states. Exclusion statistics can be 
implemented by putting n particles in m sites on a one-dimensional lattice, under 
the restriction that any two particles be at least p sites apart. For a periodic 
lattice, the number of ways of doing the above is [25] 



One can check that Di^n = •Sn(l — p) which allows to interpret the generalized 
Catalan numbers as quantum degeneracy factors. 
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Advanced Research (IFCPAR/CEFIPRA) under Project 3404-2. 




(18) 



mr(m + (1 — p)n) 



m,n — 




7 



References 



[1] G.E. Andrews, The Theory of Partitions (Cambridge University Press, Cam- 
bridge, 1998). 

[2] H.N.Y. Temperley, Statistical mechanics and the partition of numbers, the 
form of crystal surfaces, Proc. Cambridge Philos. Soc. 48, 683 (1952). 

[3] A.M. Vershik and S.V. Kerov, Asymptotics of the Plancherel measure of the 
symmetric group and the limiting shape of Young tableaux, Soviet Math. 
Dokl. 18, 527 (1977) 

[4] B. F . Logan and L. A. Shepp, A variational problem for random Young 
tableaux. Adv. Math. 26, 206 (1977). 

[5] A.M. Vershik, Statistical mechanics of combinatorial partitions and their 
hmit shapes, Funct. Anal. Appl. 30, 90 (1996). 

[6] G. Preiman, A.M. Vershik and Yu.V. Yakubovich, A local limit theorem for 
random strict partitions. Theory Probab. Appl 44, 453, (2000). 

[7] A.M. Vershik and Yu.V. Yakubovich, The limit shape and fluctuations of ran- 
dom partitions of naturals with fixed number of summands, Moscow Math. 
J. 1, 457 (2001). 

[8] A. Dembo, A.Vershik and O. Zeitouni, Large deviations for integer parti- 
tions, Markov Processes and Related Fields 6, 147 (2000). 

[9] D. Romik, Identities arising from limit shapes of constrained random parti- 
tions, preprint (2003). 

[10] L. Pastur, Lectures given at the Institut Henri Poincare trimester "Phenom- 
ena in large dimensions" (2006). 

[11] S. Shlosman, Geometric variational problems of statistical mechanics and of 
combinatorics, J. Math. Phys. 41, 1364 (2000). 

[12] P.L. Krapivsky, S. Redner and J. Tailleur, Dynamics of an unbounded inter- 
face between ordered phases, Phys. Rev. E 69, 026125 (2004). 

[13] R. Rajesh and D. Dhar, Convex lattice polygons of fixed area with perimeter- 
dependent weights, Phys. Rev. E 71, 016130 (2005). 

[14] P. Erdos and J. Lehner, The distribution of the number of summands in the 
partitions of a positive integer, Duke Math. J. 8, 335 (1951). 



8 



[15] M.E. Fisher and J.W. Essam, Some cluster size and percolation problems, 
J. Math. Phys. 2, 609 (1961). 

[16] P. Hilton and J. Pedersen, Catalan numbers, their generalization and their 
uses. Math. Intell. 13, 64 (1991). 

[17] A. Comtet, S.N. Majumdar and S. Ouvry, Integer partitions and exclusion 
statistics, J. Phys. A: Math. Theor. 40 (2007) 11255 

[18] F.D.M. Haldane, Fractional statistics in arbitrary dimensions: a generaliza- 
tion of the Pauh principle, Phys. Rev. Lett. 67, 937 (1991). 

[19] A. Dasnieres de Veigy and S. Ouvry, Equation of state of an anyon gas in 
a strong magnetic field, Phys. Rev. Lett. 72, 600 (1994); One dimensional 
statistical mechanics for identical particles-the Calogero and anyon cases. 
Mod. Phys. Lett. B 9, 271 (1995). 

[20] S.B. Isakov, Fractional statistics in one dimension. Int. J. Mod. Phys. A, 9, 
2563 (1994). 

[21] Y.S. Wu, Statistical distribution for generalized ideal gas of fractional statis- 
tics particles, Phys. Rev. Lett. 73, 922 (1994). 

[22] M.V.N. Murthy and R. Shankar, Thermodynamics of a one dimensional ideal 
gas with fractional exclusion statistics, Phys. Rev. Lett. 73, 3331 (1994 ); 
Haldane statistics and second virial coefficient ihid 72, 3629 (1994). 

[23] V. Pasquier, A lecture on the Calogero Sutherland models, in Lecture Notes 
in Physics (Integrable Models and Strings) 436, pp. 36-48 (Springer Berlin/ 
Heidelberg 1994). 

[24] A.P. Polychronakos, Physics and mathematics of Calogero particles, J.Phys. 
A 39, 12793, (2006). 

[25] A.P. Polychronakos, Probabilities and path-integral realization of exclusion 
statistics, Phys. Lett. B 365, 202 (1996). 



9 



